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We study the thermal transport properties of several quan- 
tum spin chains and ladders. We find indications for a diverg- 
ing thermal conductivity at finite temperatures for the models 
examined. The temperature at which the non-diverging pref- 
actor fc'*''-' (T) peaks is in general substantially lower than the 
temperature at which the corresponding specific heat cv(T) 
is maximal. We show that this result of the microscopic ap- 
proach leads to a substantial reduction for estimates of the 
magnetic mean-free path A extracted by analyzing recent ex- 
periments, as compared to similar analyses by phenomenolog- 
ical theories. 

PACS numbers: 



Introduction - The nature of thermal transport in 
systems with reduced dimensions is a long standing prob- 
lem and has been studied intensively in classical^ sys- 
tems either by direct numerical out-of-equilibrium simu- 
lations^ or by investigation of soliton-soliton scattering 
processes'^. There has been, on the other hand, con- 
siderably less progress for quantum systems. Huber, in 
one of the firsts works on the subject^, evaluated the 
thermal conductivity k(T') for the Heisenberg-chain with 
an equation-of-motion approximation and found a finite 
k{T), a result which is, by now, known to be wrong. It 
has been shown recently^, that the energy-current op- 
erator commutes with the hamiltonian for the spin-1/2 
Heisenberg chain. The thermal conductivity is conse- 
quently infinite for this model. The intriguing ques- 
tion, "under which circumstances does an interacting 
quantum-system show an infinite thermal conductivity" , 
is to-date completely open, the analogous question for 
classical systems being intensively studied^. 

A second motivation to study k(T) for quantum spin- 
systems comes from experiment. An anomalous large 
magnetic contribution to k has been observed^ for the 
hole-doped spin-ladder system Sri4_2.Caa;Cu2404i. For 
Ca9La5Cu2404i, which has no holes in the ladders, an 
even larger thermal conductivity has been measured^ 
raising the possibility of ballistic magnetic transport lim- 
ited only by residual spin-phonon and impurity scatter- 
ing. There is, however, up-to-date no microscopic calcu- 
lation for the thermal conductivity of spin-ladders. 

In this Letter we present a finite size analysis of the 
thermal conductance in Heisenberg Ji — J2 chains and 
ladders suggesting ballistic thermal transport in these 
two families of spin models. Following this result we 
propose a possible scenario that would account a num- 
ber of anomalies found in the thermal conductivity in 
spin systems. We apply these ideas to the case of 
Ca9La5Cu2404i, obtaining good agreement with recent 



experimental measurements. 

Models - We consider quasi one-dimensional systems 
for which the hamiltonian takes the form H = X]x=i -^x, 
where L is the number of units along the chain. We will 
consider two models, the isotropic Heisenberg chain with 
dimerized nearest and homogeneous next-nearest neigh- 
bor exchange couplings, 

Hi^'''> = j{(l + (5(-inS,.S,+i-|-aS,-S,+2} , 
and the two-leg Heisenberg-ladder: 
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77 (c/i) jj^ag been proposed to model the magnetic prop- 
erties of the spin-Peierls compound CuGeOa. The non- 
frustrated dimerized Heisenberg chain (q — 0) models 
the magnetic behavior^ of (VO)2P207. Ca9La5Cu2404i 
contains doped spin-chains and undoped spin-ladders^, 
described by H^^'"i\ 

Method - Our goal is to make connection between the 
microscopic transport-properties of the low dimensional 
models presented above and the experimental measure- 
ments of thermal conductivity mentioned in the intro- 
duction. 

The thermal conductivity is defined as the response 
of the energy-current density to a thermal gradient 
= — kVT and it has units of [k] — , where 

d is the dimension (d = 3 for experiments). The spe- 
cific form of the energy-current associated with a given 
hamiltonian is determined by the continuity equation for 
the energy density, + ' Jx = 0, which leads via 

V-Jx^= 0x^+1 -Jx^)/Ct09 



[i/x, -ffx+l] + [Hx, -ffx+2] + [-ffx-1, ^^x+l] 
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where c is the lattice constant along the chain-direction. 

In the absence of applied magnetic fields spin-inversion 
symmetry holds and the Kubo formula for k — 
lim^j^o reduces to* 



lim 
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where — J2x Jx is the total energy-current, f3 — 
1/(A;bT) is the inverse temperature and s ^ the usual 
convergence factor. 

We will examine here the possibility of infinite intrinsic 
thermal conductivity. For such a system heat transport 
is ballistic and energy is transported without dissipation. 
As an example we consider the XXZ chain with 
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a = = S and a spin-anisotropy). The total thermal 
current commutes with the hamiltonian in this case^ and 
the current- current correlation function is therefore time- 
independent. Eq. (1) reduces then to the static expecta- 
tion value k{oj = 0) = (3'^{{J^)^)/{Ls), a quantity which 
can be evaluated by Bethe-Ansatz^°. 

In general we have [J^ , H] j^O, but the thermal trans- 
port will be ballistic if 
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is finite in the thermodynamic limit L oo. The thermal 
conductivity k = k^*^^ / s then diverges for s — > 0; as it 
does for the XXZ chain. 

In reality, spin-hamiltonians like iJ^'^'*) and H^'"'^^ are 
coupled to an external environment, e.g. to phonons or 
impurities. Here we consider the case where this coupling 
is small. This coupling will then result in a finite external 
lifetime r = l/s for the cigenstates of the spin- model 
(which is assumed to be energy independent in Eq. (1)). 
Using the relation A = VsT in between the external mean- 
free path A, the spin-wave velocity Vs and r, we propose 
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to hold for the thermal conductivity k. 

Let us discuss the relation of Eq. (3) to the usual phe- 
nomenological formula^^ (here in one dimension) 



(4) 



where cy is the specific heat. When applied in order 

to analyze experimental data, the microscopic equation 
(3) will in general yield different values for the magnetic 
mean-free path A. On the other hand, we might expect 
Eq. (4) to hold at low temperature for the Hcisenbcrg- 
chain in the gapless phase, when the Luttinger-liquid 
quasiparticles are well defined and a Boltzmann approach 
is justified. Consequently we expect X/vg = cyVaX 
in this case, i.e. we expect n!-^^'' jcy = v1 in the limit 
T ^ 0. All three quantities in this equation («;(*''', cy 
and Vs) can be computed by Bethe-Ansatz for the XXZ- 
chain. One finds that this equation is exact^° in the limit 
0. 

Numerics - The computation of k^*'') (T) demands the 
whole spectrum of the hamiltonian, restricting the max- 
imum lattice size and a careful finite-size analysis is re- 
quired. In general one finds, e.g. for the XXZ-model for 
which the exact Bethe-Ansatz result is known^*^, that the 
chains with odd(cven) number of sites chains is 
a upper(lower) bound of the exact result. The value of 
k^*'^\T) decreases (increases) for chains with odd (even) 
number of sites when the lattice is enlarged. This obser- 
vation has led us to consider the average 
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of finite-size data, where Lq = 2i and Li = 2i + 1 and 
Lcff = (-^0 + Li)/2. We find that (Lefr) converges 
to the thermodynamic limit somewhat faster than tak- 
ing the limit for even (or odd) number of sites only, see 
Fig. (1). This average technique is reminiscent of the 
boundary-condition integration technique^^ for exact di- 
agonalization studies. 

The situation is similar when boundary conditions are 
changed from periodic to anti-periodic and for the ladders 
when the number of rungs changes from odd to even. We 
also find that (5) produces excellent results when applied 
to the specific heat, as illustrated in Fig. (2), where we 
compare the exact diagonalization data for ladders with 
Quantum-Montc-Carlo results^^. 

Results - After these technical remarks we discuss now 
the exact diagonalization results for the temperature de- 
pendence of The general features in all the models 
studied are: (i) A finite value of k^*'') at any finite tem- 
perature that does not vanish when the extrapolation to 
infinite lattice size is taken, (ii) A single maximum Kmax 
at an intermediate temperature Tmax{i^) smaller than the 
maximum in the specific heat Tmax{Cv)- (ni) The high 
temperature regime follows in general the law C{L)/T'^, 
as expected from Eq. (2). 

The results for are presented for J±_ = 2J|| in 

Fig. (2). We notice that the results for the specific heat 
are already converged nicely, albeit the small effective 
chain length used. We believe that the results for k^*^^ 
shown in Fig. (2) to be accurate to about 10%, enough to 
allow a detailed analysis of the experimental data which 
we will perform further below. In the inset of Fig. (2) we 
present a blowup for the data for at low tempera- 
tures for L = 3, . . . , 7. We note the systematic increase 
of K^*'') with increasing L for the odd values L = 3, 5, 7, 
indicating a finite value in the thermodynamic limit. The 
even values i = 4, 6 seem to constitute upper bounds to 

In Fig. (3) we present the results for H^'^'^\ a = 0.35 
and S — 0. We have also studied the dimcrized phase 
with S > and found similar behaviors. The maximum 
value of K^*'*) is nearly size-independent for L=8-14. The 
raise of with increasing system size L (as indicated 
by the arrows in Fig. (3)) in the low- temperature regime, 
is consistent with the notion of a finite k^*'*^ in the ther- 
modynamic limit. 

The way finite size chains approach the L ^ oo limit 
is almost identical for frustrated chains and the ex- 
actly solvable a = case, as exemplified by the size- 
dependence for the prefactor of the leading l/T^-term 
presented in the inset of Fig. (3), showing the reliability 
of the finite size analysis. We can conclude that the frus- 
tration produces a substantial drop in the extrapolated 
finite value of the prefactor. 

Analysis - Comparison with experimental results for 
the thermal conductivity ^(''^p) can be made using the 
dimensional analysis: 



(5) 
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,(e.p) = ^(^)(^)(^) k^^^) (6) 

Tic \ abc J \c J \ hvs J 

where the quantities in the brackets are dimensionless. 

is the dimensionless thermal conductance (2) which 
we will evaluate by exact diagonalization. a, b and c arc 
the lattice constants (the chains run along the c direction) 
and Nc is the number of chains per unit cell. Note that 

ic) (?ir) = ^^'^ (6) ^® ^^'^ *° extract the 

lifetime r (in units of the coupling-constant J) directly 
by comparison with the experimental k('^^p) . 

We now 
analyze the experimental data'^ for Ca9La5Cu2404i. A 
quantitative good description of the magnetic excitations 
in LagLa5Cu2404i can be obtained by F^'^'') with the 
inclusion of an additional ring-exchange term^^ (a 15% 
correction) and Jj_ ~ J||^^'^^. Here we disregard the pos- 
sible ring-exchange term and use J = J± = = 832.6 K, 
which we extracted by fitting the magnetic contribu- 
tion^ k'-'^^p^T) by 444.2/(1.0 -|- 0.093 * exp(419.6/T) + 
(T/156.4)^) (a very good fit at all temperatures, see inset 
of Fig. (4)). The gap A for the isotropic Heisenberg lad- 
der is^*' A = 0.504* J±, which leads for Ca9La5Cu2404i 
to J = 419.6/0.504 K = 832.6 K. 

Using the appropriate lattice constants and Nc = 14 
(rungs per miit cell^^) we used Eq. (5) and (6) for 
2x6 and 2x7 ladders in order to extract the spin- 
environment relaxation time r, see Fig. (4). We find 
for the lifetime t (in units of the coupling-constant 
J): Iim + i^f)- The assumption of a 

weak spin-environment coupling entering Eq. (3) is there- 
fore justified in the experimental relevant temperature 
regime. The lifetime is, to give an example, 132 times 
larger than the coupling constant J at T = 150 K, lead- 
ing to an external broading of the energy levels of only 
1/132 = 0.76%, in units of J. 

To estimate the effective mean-free path we use 
A(r) = Vs{T)t(T), were we have used for Vg the ther- 
mal expectation value of the magnon- velocity within the 
independent-triplet model: 

If 1 
= — / dk e' {k) n{e{k)) = — / dsn{e) 
Z J ^ Ji\ 

where Z = J dkn{e{k)) is the partition function and 
n{e{k)) = 3/(exp(/3£(fc) +3). We approximated the one- 
magnon dispersion e{k) by 2e^{k) = (A^ -|- A|) + (A| — 
A2)cos(fc). For Jj_ = J|| we have^^ that A2 w 3.8A. 
In the low-temperature limit VsiT) ~ VT holds and the 
mean-free path A(T) is therefore less divergent at low- 
T than the relaxation time t(T) see Fig. (4). At very 
low temperature we expect impurity-scattering to be- 
come relevant and A(T) to plateau-off. 

Discussion - Our results for the mean-free length A(T) 
for Ca9La5Cu2404i obtained by the microscopic formula 
(3) are substantially smaller than the one obtained us- 
ing (4). At T = 100 K Hess et af estimated a large 



A « 3OOOI, in contrast to our result A(IOO) « 176A. 
This is, in a certain sense, surprising, since the effective 
A obtained from (4) might be thought to contain addi- 
tional spin-spin scattering. Physically, the reason for our 
smaller mean- free path stems from the fact that k^*^^ (T) 
peaks at substantially smaller temperatures than cy(T), 
see Fig. (3). We note, that all excitations contribute 
to cy on an equal footing, but that magnetic excita- 
tions near the bottom of the one-magnon band seem to 
contribute dominantly to the thermal conductivity, lead- 
ing to a substantial reduction of the temperature where 
i^i^^){T) is maximal, with respect to cy(T). 

Conclusions - We have presented a microscopic ap- 
proach of thermal transport in quasi one-dimensional 
spin models. For several models we find indications for a 
diverging thermal conductivity. A finite thermal conduc- 
tivity is obtained when couplings to external degrees of 
freedom are taken into account with the relaxation-time 
approximation. We have analyzed recent experiments 
for spin-ladder compounds and found only weak spin- 
environment coupling for the ladder compounds, which 
decreases fast with decreasing temperature. Our esti- 
mates for the mean-free length turned out to be sub- 
stantially smaller than previous estimates using simple 
phenomenological formulas. These results highlights the 
importance of microscopic theories for transport in quasi 
one-dimensional quantum-spin systems. 
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FIG. 1. k'*''', as defined by Eq. (2), for the Heisen- 
berg-chain as a function of temperature. Inset: Data for 
L = 5 ... 14, in comparison with the Bethe-Ansatz re- 
sult^" (full line). Main panel: Results using Eq. (5), with 
Le// = (-£-0 + Li)/2. We note, that the position of the maxi- 
mum in k'*'''' can be determined confidently with the averag- 
ing procedure. 
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FIG. 2. Temperature dependence of and the dimen- 

sionless specific heat for ladders with twisted boundary con- 
ditions and J_L = 2J||. Both magnitudes are computed using 
the average procedure (5) described in the text. The specific 
heat for a 2 x 100-ladder computed with QMC is plotted for 
comparison (squares). Inset k'"'' for L=3,5,7 (from bottom 
up) and L = 4, 6 (from top down) in the low-temperature 
region. 




FIG. 3. k'-*''-' and dimensionless specific heat for the frus- 
trated Heisenberg chain with a = 0.35 and different lattice 
sizes L=8,10,12,14, the arrows indicate increasing L. Inset : 
Finite size analysis (L=6...14) of the high temperature residue 
C{L) defined as k'*'*' = C{L)/T'^ at T > J for different val- 
ues of a. 
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FIG. 4. t(T) (right-scale) and A(r) (left-scale) for 
CagLa5Cu2404i extracted using Eq. (6). The experimental 
data by Hess et al7 is shown in the inset (filled circles) to- 
gether with an analytic fit (solid line) discussed in the text. 
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